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ABSTRACT 


A three-dimensional Navier-Stokes code has been de- 
veloped for analysis of turbomachinery blade rows and 
other internal flows. The Navier-Stokes equations are 
written in a Cartesian coordinate system rotating about 
the i-axis, and then mapped to a general body-fitted 
coordinate system. Streamwise viscous terms are ne- 
glected using the thin-layer assumption, and turbulence 
effects are modelled using the Daldwin-Lomax turbulence 
model. The equations are discretized using finite differ- 
ences on stacked C-type grids and are solved using a mul- 
tistage Runge-Kutta algorithm with a spatially-varying 
time step and implicit residual smoothing. 

Calculations have been made of a horseshoe vortex 
formed in front of a flat plate with a round leading edge 
standing in a turbulent endwall boundary layer. Compar- 
isons are made with experimental data taken by Eckerle 
and Langston for a circular cylinder under similar con- 
ditions. Computed and measured results are compared 
in terms of endwall flow visualization pictures and total 
pressure loss contours and vector plots on the symme- 
try plane. Calculated details of the primary vortex show 
excellent agreement with the experimental data. The cal- 
culations also show a small secondary vortex that was not 
seen experimentally. The calculations required about 1.6 
million words of storage and 1.3 hours of CPU time on a 
Cray X-MP computer. 

Calculations have also been made of an annular tur- 
bine stator that has been tested experimentally at NASA 
Lewis. The Mach number ranged from from about 0.21 at 
the inlet to 0.67 at the exit, and the axial chord Reynolds 
number was about 1.7 x 10 5 . Computed surface pressure 
distributions compare well with measured values at three 
span-wise locations. The endwall boundary layers pro- 
duce horseshoe vortices at the leading edge of the blade. 
Computed wake profiles resemble the measured profiles, 
but computed efficiencies are lower than measured values 
by a factor of two. 


INTRODUCTION 

Much of our work in the past few years has in- 
volved the analysis of two-dimensional or quasi-three- 
dimensional blade-to-blade flows in turbomachinery 1,2 . 
In these two references, both Euler and Navier-Stokes re- 
sults were presented for each blade row considered, and in 
each case a significant viscous effect was observed. These 
effects included such things as a pressure-surface or ‘cove’ 
separation on an axial turbine blade, a reduction in the 
choking mass flow and a change in shock location for an 
axial compressor blade, and a reduction in the peak Mach 
number and shock strength near the leading edge of a 


centrifugal impeller. Each of these two-dimensional vis- 
cous effects can be expected to lead to secondary flows in 
three dimensions. Complex geometries, endwall bound- 
ary layers, tip clearance effects, etc. also lead to three- 
dimensional flows in turbomachinery. It is the goal of 
the present work to begin to predict some of these three- 
dimensional viscous effects. 

Several steady three-dimensional analyses for tur- 
bomachinery have been published lately. Among them 
are the work of Dawes 3 and Subramanian and Bozzola 4 . 
Both used Runge-Kutta schemes implemented on sheared 
H-type grids. H-grids are particularly easy to generate 
and implement for turbomachinery calculations, but suf- 
fer from poor leading-edge resolution. Rai 5 has published 
a notable analysis of unsteady three-dimensional rotor- 
stator interaction in an axial turbine. He used a third- 
order accurate upwind implicit scheme and a system of 
patched and overlaid O- and II-type grids for good reso- 
lution of viscous phenomena. His analysis has not been 
used for steady flows, however. 

In this paper we describe a numerical method for 
analyzing three-dimensional viscous flows in isolated tur- 
bomachinery blade passages. The underlying Cartesian 
formulation allows the method to be applied easily to 
both Cartesian and cylindrical geometries. Stacked C- 
type grids give good resolution of critical leading-edge 
regions. 

A multistage Runge-Kutta scheme is used to solve 
the finite-difference form of the thin-layer Navier-Stokes 
equations with a Baldwin-Lomax turbulence model. A 
spatially-varying time step and implicit residual smooth- 
ing are used to accelerate convergence of the scheme to a 
steady state. Two calculations are presented to validate 
the analysis. The first calculation shows the formation 
of a horseshoe vortex at the leading edge of a flat plate 
with a round leading edge that stands in a turbulent end- 
wall boundary layer. Comparisons made between calcu- 
lated results and experimental data for a circular cylin- 
der under similar conditions show excellent agreement 
between static pressure distributions and flow visualiza- 
tion pictures on the endwall, and between static and total 
pressure contours and velocity vector plots made on the 
symmetry plane. The second calculation is of the flow 
through an annular turbine stator. Comparisons made 
between calculated and measured static pressure distri- 
butions compare well at three span-wise locations. 


GOVERNING EQUATIONS 

The Navier-Stokes equations are written in a Carte- 
sian (x, y, z) coordinate system rotating with angular 
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velocity Cl about the x-axis. The rotation introduces 
source terms in the y- and z- momentum equations. The 
Cartesian equations are mapped to a general body-fitted 
(£, 'll () coordinate system using standard techniques. 
The ^-coordinate direction is assumed to follow the flow 
direction and the thin-layer approximation is used to 
drop all viscous derivatives in this direction. All vis- 
cous terms in the cross-channel (rj, f) plane are retained, 
with an option in the code to delete all cross-derivatives 
if desired. The resulting equations are as follows. 

d%q -I- J\d^P 4- df)P 4- — Re ^{dijPv 4- d f (3V)] = H 

(!) 

where: 

q=J~ 1 \p, pu, pv, pw, ef 
= 1 0, 0, -Clpw, Clpv, Of 

= J~ l X 

\pU', puU' 4- fxP, pvW 4- £ v p, pwU‘ 4- £,p, eU‘ 4- pU} T 

P = J _1 x 

| pV', puV 1 4- T} x p, pvV' 4- r] v p, pwV' 4- f UP, eV‘ + pV f 
G= J~ l x 

| pW\ puW' 4- < z p, pvW' + fy p, pwW + ( x p, eW‘ + pW f 

( 2 ) 

The velocities in q are absolute with respect to the co- 
ordinate system fixed to the blade. Relative velocities 
(denoted by a prime) are given by: 

u' = u 
v' = u — Clz 
w 1 = w 4- Cly 

( 3 ) 

and the relative contravariant velocitiy components are 
given by: 


U' = £xU 4- Zy v ' 4- £ t w' 

V‘ = r) z n + rj y v‘ + r\ x w‘ 

W' = fxU 4- C u v‘ 4- $ x w' 

Note that although u' = u, U' ^ U . 

The energy and static pressure are given by: 


( 4 ) 


e = p [C V T + (u 2 4- v 2 4- u: 2 )/2] (5) 

P = (-7 — 1) [« — P(“ 2 4- u 2 4- tu 2 )/2] (6) 


and 


Ci = r, 2 + rj 2 + r, 2 

Ci = |(f?x9 0 u 4- ri u d n v 4- »J*dqU)) 

C 3 = rjxfx + n v Cv + »7.fx 

Ci = \(c x d < u 4- fsAu 4- CtdfW) 

Cs = rj x d f u + ri v d ( v + r) x d ( w 

(7c) 

Terms multiplied by C x and C 2 lead to non-mixed second 
derivative viscous terms like u nn , while terms multiplied 
by C 3 -C 3 lead to mixed-derivative terms like u nf . The 
viscous flux vector G can be written similarly, with di- 
rections rj and f everywhere interchanged. 

Metric terms are defined using the following rela- 
tions. 


Zz fix C* 

Zy Vy Cy 

. £» ?* . 


= J 


Vrt z < - y f z„ 
x t z n ~ x n z r 
x r,y< - x (Vn 


where 


V< z i ~ Vi * f 
X ( Z ! ~ x { z ( 
* f y« - x iV t 


yt z n - Vn z i 

x r) z t ~ x i z n 

x (y n - x r,yt . 


(8) 


j = 

(x t y n z f 4- x f y £ z„ 4- x n y ( z t - x £ y f z„ - x n y £ z f - x f y„z f ) -1 

( 9 ) 

The equations are nondimensionalized by arbitrary 
reference quantities (here the inlet total density Poref 
and the total sonic velocity c 0re / were used,) and the 
Reynolds number Re and Prandtl number Pr are defined 
in terms of these quantities. The equations assume that 
the specific heats C p and C v and Prandtl number are 
constant, that Stoke’s hypothesis is valid, and that the 
effective viscosity for turbulent flows may be written as 


Me// — Mlam "h Mtur6 (10) 

where the laminar viscosity is calculated using a power 
law function of temperature: 


plam = ( T \ n 
Mre/ \ T re / / 


(U) 


with n = | for air. 


TURBULENCE MODEL 


Using Stoke’s hypothesis, A = — |p, the viscous flux 
Pv can be written as follows: 


P v = F 2 , F 3 , F t , F & ] r (7a) 


where 

F 2 = C x d n ii 4- Cjfh 4- C 3 d ( u — Cm, 4- C 3 ( x 
F 3 = C x d n v 4- Cirjy 4- C$d ( v — CiTjy 4- Csf y 
Fi = C\d n w 4- Citfz 4- CsdfW — Cirj x 4- C&£ c 

F s = ^(Cid^CyT) 4- C 3 d ( (C v T)) + uF 2 + vF 3 + wF 4 

(76) 


The Baldwin-Lomax algebraic two-layer eddy viscos- 
ity model 6 is applied on cross-channel (q, f) planes. Two 
modifications to the standard model are made to account 
for the endwall boundary layer and the blade boundary 
layer and wake, and their interaction in corners. 

First, the distance from the wall is calculated using 
the Buleev 7 length scale d: 


d = 


2 s n S( 


*n + *( + (* 2 + s ?) 1/2 


( 12 ) 


where s n and s ( are normal distances from the walls in 
the tj- and f- directions respectively. This length scale 
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has the desirable property that d approaches the normal 
distance from one wall at large distances from the other 
wall. 

Secondly the turbulent viscosities are calculated 
across each boundary layer or wake separately, then the 
total turbulent viscosity is taken as the vector sum of the 
components, Le., 


Pturb i^turb rj Pturb (•) (^) 

This assumption has the desirable properties that outside 
of one viscous layer Pturb takes on values calculated for 
the other layer, that it goes to sero in the core flow, and 
that near corners it accounts for both walls. 


COMPUTATIONAL GRID 


For subsonic outflow the exit static pressure is spec- 
ified and (p, pu, pv, pui) are extrapolated. For Cartesian 
geometries the exit pressure is constant. For annular 
geometries the hub pressure is specified and the radial 
pressure distribution is found by integrating the axisym- 
metric radial momentum equation: 


dp 

dr 



(17) 


Sidewalls and the trailing-edge cut are treated as 
periodic boundaries. 

On the blade surface V' = 0, and for viscous flows 
U' = W' = 0. Blade surface pressures are found from the 
normal momentum equation. On the hub (f = 1) and tip 

(f = fmax)* 


(?x£x + Svtv + (fx'Jx + fj/fJw + 


Two-dimensional body-fitted grids for this work 
were generated using the GRAPE code developed by 
Sorenson ®. Three-dimensional grids were formed by 
stacking the 2-D grids. Figure 1 shows a 3-D grid around 
a plate with a round leading edge. For annular geome- 
tries the 2-D grids were stacked along a radial stacking 
line and stretched in the rj-direction so that the blade 
shape remained constant and the angular pitch of the 
outer (periodic) boundary remained constant. 


+(f* + f* + R p = -p[ n(fj,w - ? r v) 

+U'(Szd(U+<v d e' 1 +$‘ a ii w ) + V'(Sxd n u+$ u d n v+c e 9 n w)\ 

(18) 

On the blades ( rj = 1) the normal momentum equation 
can be found from (18) by replacing f everywhere by rj 
and V' by W 1 . 


MULTISTAGE RUNGE-KUTTA ALGORITHM 


BOUNDARY CONDITIONS 

At the inlet, total temperature T 0 rc j is specified as 
a constant. A f -distribution of total pressure (Po/Pq rc /) 
is specified, c/2/2- as a constant or as appropriate for an 
inlet boundary layer with given thickness and a power- 
law velocity profile. For Cartesian geometries the (z, y) 
and (z, z) flow angles are specified. For cylindrical ge- 
ometries the (z, y) flow angle is replaced by the inlet whirl 
rvg. 

For supersonic inlet flows, all flow variables are spec- 
ified at the inlet. For subsonic flows the inlet condi- 
tions are updated each iteration by extrapolating the 
upstream-running Riemann invariant R~ based on the 
absolute total velocity Q = \Ju 2 + v 2 + w 2 to the inlet. 

R~ - Q —r (14) 

7—1 

The total velocity is then found from the To re/ using: 

( 7 - 1 )R~ + v/2(l - 7)(*-) 2 + 4(7 + 1 )C P T 0 rcf 

V fr+1) 

(15) 

Velocity components are found from Q and the specified 
angles or whirl. Within the endwall boundary layer, that 
is, where Po/Porc / < li the v and w velocity components 
are found by extrapolation from upstream. The density 
is found using: 

P _ Po 

PO re/ Po ref 


The governing equations are discretized using a 
node-centered finite difference scheme. Second order cen- 
tral differences are used throughout. 

The multistage Runge-Kutta scheme developed by 
Jameson, Schmidt, and Turkel 8 is used to advance the 
flow equations in time from an initial guess to a steady 
state. If we rewrite (1) as 

dtq^-JlRt-iRv + D)} (20) 

where Rj is the inviscid residual including the source 
term, Ry is the viscous residual, an D is an artificial 
dissipation term described in the next section, then the 
multistage Runge-Kutta algorithm can be written as fol- 
lows: 

90 = <7n 

91 = 9o — aiJAt\Ri q 0 - (Rv + D) 9o] 

9* = <?o - a k JAt\Ri q k - j - {R v + D) q 0 ] 

9n+i = 9 k 

For efficiency both the physicial and artificial dissipation 
terms are calculated only at the first stage, then are held 
constant for subsequent stages. 

ARTIFICIAL DISSIPATION 

The dissipative term D in (20) is a nonconservative 
version of that used by Jameson et al. 9 It is given by: 

Dq = (Z>£ + D n + D ( ) q 


7-1 (Q-R-) 


(16) 
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c 0 ref \ 


(22a) 



where the £ -direction operator is given by 
D ( q = C(V 2 q ti - V t q u(t ) 


Jameson 11 . The technique involves replacing the residual 
calculated in (20) with a value that has been smoothed 
(226) by an implicit filter, Le., 


where 



(22c) 


is a coefficient that cancells similar terms in (21). To 
minimise the artificial dissipation in viscous regions we 
reduce C linearly across several grid points to sero at the 
walls. 


The terms Vj and V* are given by: 


V 2 = ntmax(ui +1 ,v i ,u i - 1 ) 

V< = max (0, /i 4 - V 2 ) 


where 


and 


_ Ifi+i.y ~ 2 -P«.y + 

1,1 |Pi+i,/ + 2P<,y + P,_i,y| 


(22d) 

(22c) 


M2 = 

A»« = 

( 22 /) 

In smooth regions of the flow the dissipative terms are of 
third order and do not detract from the formal second- 
order accuracy of the scheme. Near shocks i/ ti y is large 
and the dissipative terms become locally of first order. 



THREE-DIMENSIONAL STABILITY LIMIT 

Applying a linear stability analysis to the inviscid 
form of (20-21) gives the following expression for the time 
step. 


At < 


where 


CFL 

1.M + f„M + i.M + 

(23 a) 


l x = |£*| + In* | + k* | 

L = IGI + ln»| + Ifvl 

'* = |(*| + 1»7*| + k*| 

(236) 

The Courant limit for a particular multistage scheme de- 
pends on the number of stages and the choice of coeffi- 
cients cti . See Ref. 9 for several examples. 

To accelerate convergence to a steady state we use 
the maximum permissible time step at each grid point so 
that the Courant number is constant everywhere. The 
time step is calculated once based on the initial condi- 
tions. It is stored and is not updated during the calcula- 
tions. 


IMPLICIT RESIDUAL SMOOTHING 

Residual smoothing was introduced by Lerat (see 
for example Ref. 10) for use with the Lax-Wendroff 
scheme and was later applied to Runge-Kutta schemes by 


(1 - e e « £e )(l - £„«„)( 1 - e t 6")R = R (24) 

where 6 nn , and S (( are standard second difference 
operators and £(, e n , and e f are smoothing parameters. 

Linear stability analysis has shown that the Runge- 
Kutta scheme with implicit residual smoothing may be 
made unconditionally stable if the e smoothing parame- 
ters are made sufficiently large. In one dimension 



gives unconditional stability if A* is the Courant limit 
of the unsmoothed scheme, and A is a larger operating 
Courant number. In three dimensions different e’s may 
be used in each direction, and their magnitudes may be 
often be reduced below the value given by Eq. (25.) 


RESULTS 

Two sets of computed results are presented for pre- 
liminary validation of the code described above. The first 
set of results shows the structure of a horseshoe vortex 
formed at the base of a cylinder standing in a turbulent 
boundary layer. The second set of results is for turbu- 
lent flow through an annular turbine cascade. Computed 
results are compared to experimental data in each case. 

When a boundary layer approaches a local obstruc- 
tion such as the leading edge of a turbine blade, the low 
momentum fluid in the boundary layer often cannot over- 
come the local pressure gradient and the flow separates 
from the wall. In front of the obstacle the separation 
creates a vortex which convects around the sides of the 
obstacle and forms a characteristic horseshoe-shaped flow 
region. 

Eckerle and Langston 12 have made detailed mea- 
surements of the horseshoe vortex in front of and around 
a cylinder of diameter D centered between the sidewalls 
of a wind tunnel. Test conditions included an inlet Mach 
number of 0.084, Rep = 5.5 x 10 5 , and an upstream tur- 
bulent boundary layer thickness 6 = 0.1D. Detailed sur- 
face flow visualisation and static pressure measurements 
were made, and static and total pressure measurements 
were taken using a five-hole probe. 

Figure 1 shows the grid used to compute Eckerle 
and Langston’s flow. The grid shown has been coars- 
ened for clarity. The actual grid had 65 x 49 x 25 points 
with an initial spacing at the walls As,- == 0.001D. To 
avoid questions of trailing-edge vortex shedding, a tail 
board was added from the back half of the cylinder to 
the exit boundary. The base grid is approximately 6D 
square and .5D high, to match the dimensions of Eck- 
erle and Langston’s wind tunnel test section. We used a 
symmetry condition at mid span but computed the full 
symmetric flow side-to-side. 

The experimental inlet Mach number of 0.084 is too 
low for the compressible algorithm used here, so the cal- 
culations were run with M,„ = 0.2. The peak Mach 
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number was 0.34 at the cylinder-plate junction, so the 
flow was essentially incompressible. The Reynolds num- 
ber and inlet boundary layer thickness were matched to 
the experimental data. 

Even with M; n = 0.2 the four-stage Runge-Kutta 
scheme seemed to converge poorly, so we eventually ran 
the calculations with a two-stage scheme with a,- = 
(1.2,1.), CFL = 4., and implicit residual smoothing at 
each stage with e ( = 2., e n = e ( = 4. The initial conver- 
gence rate was fast, as shown in Fig. 2 by the histories 
of the maximum and rms residuals, but after 1000 iter- 
ations the solution showed only a tiny horseshoe vortex 
that did not match the data. Over the next 3000 iter- 
ations the residuals changed little, but the vortex grew 
and moved upstream until it stabilized at the position 
shown later. The total solution took 1.3 hours on a Cray 
X-MP computer. 

Figures 3 and 4 compare the calculated and experi- 
mental data on the endwall. Figure 3 shows contours of 
constant static pressure coefficient, defined by: 




P. ~ P. ref 

Po re/ — Pi ref 


(26) 


Computed contours are on the bottom and measured 
contours are on the top. The calculations show excel- 
lent agreement with the data from the symmetry plane 
to about 45 degrees around the cylinder, where the in- 
fluence of the tailboard becomes apparent. The lower 
part of Figure 4 shows a flow visualization picture made 
with ink dots on the endwall. It clearly shows the sepa- 
ration line and reverse flow region ahead of the cylinder. 
The computed vector plot at the top of the figure shows 
close agreement with the measured separation line loca- 
tion and flow directions. The vectors are one point off 
the endwall, and are all drawn to the same length, so 
they show direction only. 

Figure 5 compares the static pressure coefficient dis- 
tributions on the endwall along the symmetry line ahead 
of the cylinder. The experimental data (circles) show 
a general pressure rise upstream due to the cylinder 
blockage, but also a large dip in the pressure inside the 
separated region. Two vertical bars indicate substan- 
tial unsteadiness in the experimental data. Eckerle and 
Langston included a 2-D potential solution (solid line) 
for comparison. The computed solution (line with trian- 
gles) shows good agreement with the data ahead of the 
saddle point, but falls short of predicting the magnitude 
of the pressure dip. The discrepancy may be due to the 
unsteadiness in the real flow or to lack of resolution in 
the computed solution. 

Figure 6 compares computed (bottom) and experi- 
mental (top) velocity vectors and total pressure loss co- 
efficient contours on the symmetry plane upstream of the 
cylinder. The total pressure loss coefficient is defined by: 


C pt - 


Po ref ~ Pp 
Po ref ~ P$ ref 


(27) 


These contours show nearly horizontal boundary layer- 
like flow upstream that rolls up into a vortex with a high 
loss core. The low loss fluid above the boundary layer 
curves down the face of the cylinder and carries high 
momentum fluid to the region near the cylinder-endwall 


junction. There is excellent agreement between the com- 
puted and experimental data. 

The velocity vectors show how the flow rolls up to 
form a horseshoe ahead of the cylinder. Experimental 
velocity vectors are missing in areas where the flow an- 
gle exceeded the calibrated range of the five-hole probe. 
From this data Eckerle and Langston 12 concluded that 
“The reverse flow did not roll up to form a vortex, 
however. The vectors clearly show that a closed vor- 
tex was not present in the plane of symmetry, though 
positive pitch angles in a portion of the reverse flow at 
R/D = 0.72 and Y/D = 0.02 may indicate the start of 
vortex formation. Rather than rolling up, flow passed 
out of the plane and proceeded tangentially around the 
cylinder.” 

The computed vectors clearly show a vortex in the 
symmetry plane. A small counterrotating secondary vor- 
tex is also shown at the cylinder-endwall junction. The 
dimension of the secondary vortex is about 1.5 times the 
diameter of the five-hole probe, and would have been 
nearly impossible to detect experimentally. 

The second set of results is for an annular cascade 
of constant profile turbine stator vanes developed and 
tested at NASA Lewis 13,14 . The annular ring has 36 
vanes with a hub-tip radius ratio of 0.85 and a tip di- 
ameter of 508 mm. The vanes themselves are 38.10 mm 
high and have an axial chord of 38.23 mm. Design flow 
conditions are for a fully axial inflow with a hub-static to 
inlet-total pressure ratio of 0.6705. These conditions cor- 
respond to average inlet and exit Mach numbers of 0.211 
and 0.665 respectively. The Reynolds number based on 
axial chord is 1.73 x 10®. 

A grid consisting of 97 x 31 x 33 points with an initial 
spacing at the wall of 0.0002 of a blade chord was used 
for the flow calculations and is shown in Figure 7. 

The calculation was run with a four-stage scheme 
with ct% = (1/4, 1/3, 1/2, 1) and CFL = 5.5, using 
implicit residual smoothing after each stage with = 
( n = ( ( = 0.75. Convergence histories for the annular 
cascade calculation are shown in Fig. 8 where the log of 
the maximum and rms-averaged residuals have dropped 
approximately 3.5 orders of magnitude in 1500 iterations. 
The total CPU time was approximately 54 minutes for 
this calculation. 

Mach number contours at mid span are shown in Fig. 
9 to illustrate the blade boundary layer and wake thick- 
nesses. There are approximately 12 grid points across 
the pressure surface boundary layer. 

The inlet boundary later thicknesses were specified 
as 1.9 percent span on the hub and 7.1 percent span on 
the tip, corresponding to the measured data in Ref. 13. 
In Fig. 10, velocity vectors with superimposed contours 
of total pressure show how these boundary layers roll up 
into horseshoe vortices at the leading edge of the blade. 
The primary vortices are considerably smaller than the 
inlet boundary layers, and each primary vortex has an 
even smaller counterrotating secondary vortex associated 
with it. 

Figure 11 shows a comparison between the calcu- 
lated surface static pressure distribution and data ob- 
tained by Goldman and Seasholtz 13 at locations of 13.3, 
50.0, and 86.7 percent span. The calculated results agree 
very well with the experimental data at all spanwise lo- 
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cations. Although the blade section is constant from hub 
to tip, the increased pitch at the tip increases the loading 
of the tip section considerably. 

Figure 12 compares calculated (top) and measured 
(bottom, Ref. 14) efficiency contours on a cross-channel 
surface located 1/3 axial chord downstream of the trail- 
ing edge. The kinetic energy efficiency is theoretically 
independent of the axial position and is defined by the 
following: 

V Q actual t 'Q ideal 

where 

Qideal = 2C p To (l - (28) 

and To is taken as constant. 

The efficiency contours clearly delineate the endwall 
boundary layers and wake. The computed wake is thin- 
ner than the measured wake and has higher losses (lower 
efficiency) at the center. Integrating the efficiencies over 
the entire area gives a total efficiency of 0.960 for the real 
machine and 0.923 for the calculations. The discrepancy 
could be due to inadequate resolution of the thick, round 
trailing edge, an inadequate turbulence model, or possi- 
bly to unsteady vortex shedding in the real flow. The 
high calculated losses seem not to be due to numerical 
dissipation. Indeed, the computed wake appears to be 
less dissipative than the measured wake, and the com- 
puted loss is remarkably insensitive to the artificial vis- 
cosity coefficient. Considerable work is clearly needed in 
the area of modelling the flows leaving the blunt trailing 
edges commonly found on turbomachinery blades. 


CONCLUDING REMARKS 

A numerical analysis has been developed for three- 
dimensional viscous internal flows. The analysis solves 
the 3-D Navier-Stokes equations written in a general 
body-fitted coordinate system, including rotation about 
the z-axis. The thin-layer approximation is made in the 
streamwise direction but all viscous terms are included 
in the cross-planes. The Baldwin-Lomax eddy-viscosity 
model is used for turbulent flows. 

An explicit multistage Runge-Kutta scheme is used 
to solve the finite-difference form of the flow equations. 
A variable time step and implicit residual smoothing are 
used to accelerate the convergence of the scheme. Con- 
vergence rates are slow at low Mach numbers but reason- 
able for typical turbomachinery applications. We hope 
to improve convergence rates by adding multigrid to the 
code. 

Results showing the development of a horseshoe vor- 
tex in a turbulent boundary layer ahead of a cylinder are 
presented to validate the analysis. Excellent agreement 
is found between the computed results and experimental 
data for this case. 

Results for an annular turbine cascade also show 
horseshoe vortex development. Very good agreement was 
found between measured and computed surface pressure 
distributions. Computed losses were high and show the 
need for improved modelling of flows around blunt trail- 
ing edges commonly found in turbomachinery. 


It is felt that with additional work and experience 
this analysis will prove to be a useful tool for investigat- 
ing three-dimensional viscous flow phenomena in turbo- 
machinery. 
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Figure 1. Computational grid for the horseshoe vortex 
problem. 



Figure 3. Comparison of measured (top) and computed 
(bottom) static pressure contours on the endwall. 
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Figure 5. Comparison of measured and computed static 
pressure distributions on the endwall symmetry line. 



Figure 2. Convergence history for the horseshoe vortex 
problem. 



Figure 4. Comparison of endwall flow visualisation (bot- 
tom) and computed flow direction vectors (top) ahead of 
the cylinder. 



Figure 6. Comparison of measured (top) and computed 
(bottom) velocity vectors and total pressure loss contours 
on the symmetry plane. 
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Figure 7. Computational grid for the annular turbine 
cascade. 



Figure 9. Mach number contours at mid-span. 



Figure 8. Convergence history for the annular turbine 
cascade problem. 




Figure 11. Comparison of measured and computed pres- 
sure distributions at three span-wise locations. 



Figure 12. Comparison of measured (bottom) and com- 
puted (top) efficiency contours 1/3 axial chord down- 
stream of the trailing edge. 
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